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Abstract. Among the still unmeasured neutrino properties, the third neutrino mixing angle, 
^13, is likely to be the one we will next find out. In this contribution, first a brief summary of 
the limits and the preliminary measurements of this angle is given. Second a critical assessment 
of a widely used formula connecting two- and three-fiavor evolution is provided. 



1. Introduction 

Solar, accelerator, and reactor experiments performed during the last two decades firmly 
established that neutrinos mix. In general, for three flavors, the neutrino mixing matrix contains 
three mixing angles, one CP-violating phase and two Majorana phases. The Majorana phases are 
accessible only in double beta decay experiments. Two of the mixing angles are well measured. 
Prom solar and reactor experiments we know the value of one of them (usually denoted ^12)- 
In the atmospheric neutrino experiments we measured the other one (usually denoted ^23)- 
Contrary to the original (and naive) expectations, both angles are large. In fact one gets 
^23 ~ vr/4. To have the complete information on neutrino mixing, we still need to know three 
more pieces of information: i) the value of the remaining mixing angle, ^13, ii) the mass hierarchy; 
and iii) the value of the CP-violating phase. Of these three, the remaining mixing angle is likely 
to be the next piece to be measured. 

One popular parametrization of the neutrino mixing matrix is 
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T23T13T12 = C723 ^23 1 -5i2 C12 (1) 



where Cij = cos Oij, Sij = sin 9ij, and 6c p is the CP-violating phase. Note that although the 
CP-violating phase multiplies 5i3 in Eq. ([T]), there is nothing special about this choice; one 
could have placed 5 elsewhere in the mixing matrix. In fact the combination 

P{iy^ ^ Ue) - ^ ^e), (2) 

which quantifies the CP-violation in neutrino oscillations, is proportional to the sines of all the 
mixing angles. However, since we know the other mixing angles are non-zero, observation of 
CP-violation in neutrino oscillations hinges on a non-zero value of ^13. 

The current limit on ^13 quoted by the Particle Data Group is sin^(2^i3) < 0.15 at 90% 
C.L. [T]. There are strong hints for a non-zero value of ^13. Already several years ago, it was 
pointed out that the slight discrepancy between solar neutrino experiments and the KamLAND 



reactor neutrino experiment could be explained by a small, but non-zero value of ^13 [2]. Indeed 
global fits indicate a slight preference for a non-zero value of ^13 [3] . (A recent brief summary of 
the 013 phenomenology is given in [3]). Recent low-threshold analysis of the Sudbury Neutrino 
Observatory yields ^13 = 7.2^2 8 degrees, albeit with non-Gaussian errors [5J. A global analysis 
by the KamLAND cohaboration indicates sin^ 6^13 = 0.009^0:0^7 at the 79% C.L. [6]. Currently 
three reactor neutrino experiments (Double Chooz [3, Daya Bay [8], and RENO [9j) as well as 
the T2K accelerator experiment [TO] are working on a measurement of this angle. So far the best 
hint for a non-zero value of ^13 came from the T2K experiment which announced a 2.5cr result: 
0.03 < sin^ 2^13 < 0.28 for the normal and 0.04 < sin^ 2^13 < 0.34 for the inverted hierarchies, 
both for (5cp = [llj. 



2. Reduction of three-flavor evolution 

To extract ^13 from solar neutrino data many times the formula 

Psxsit^e ^e) = COS"^ 6*13 ^2x2(1^6 Z^e with iVg COS^ 6*13) + sin'^ 6*13 (3) 

is used. Since many fits to solar neutrino data that give a non-zero value of of 6*13 use Eq. (I3|), 
here we wish to give an assessment of the accuracy of this formula. To do so, following Ref. 
|12j . we introduce the combinations 

§^ = COS023*/. -Sin023^r, (4) 
= sin023^'M + COS023*r, (5) 

which leads to the MSW evolution equation 
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H = T13T12 I ^2 Tl2Tl3 +10 0). (7) 
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In writing Eq. ([7]), since it is proportional to the identity, the neutral current contribution is 
dropped and 

Ve^,{x) = V2GfN,{x). (8) 
To write the evolution operator associated with Eq. ([6]), 

i^U = HU, (9) 
following the procedure of Ref. [13] we introduce the evolution operator U: 

U = Ti3UTl3 (10) 

which satisfies the equation 

4u = «U (11) 



with 

/iy - A2icos26li2 I A21 sin 26*12 ^Vf^eSmlOn 

T|3HTi3 = H= iA2isin20i2 -i^ + A21 cos20i2 

V lV^esm29n ^Asi + A32) + V^.e - '^V ^ 

(12) 

where we took the CP-violating phase to be zero and introduced the modified matter potential 

V = V^eCOS^0is (13) 

as well as the quantity 

mf — m3 6m% , ^ 

Ai,- = — ^ = ^. (14 



Writing the Hamiltonian in Eq. (|12p as 

U = Uq + Ux (15) 



with 



and 



iy-Aai cos 2012 |A2i sin 2012 \ 

^0=1 ^Asi sin 2012 -^F + A21 cos 20i2 , (16) 

i(A31 + A32) + F^e-iV^/ 

1 /O r 

Ux = -sin20i3l^^e I I , (17) 
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we see that Eq. (fTTj) can be solved perturbatively, assuming 0i3 is small. (In seeking such 
solutions we will take all the evolution operators at the neutrino production (t=0) to be the 
identity operator). The contributions which are lowest order in 0i3 will come from from the 
quantity V and from the transformation in Eq. (jlOp . Eq. ()16p implies that if 0i3 were zero, ^'t- 
would decouple from the other flavors. Writing U = UqUi, the solution of 

i^Uo = ^oUo (18) 

is 

where S is the evolution operator for the two-flavor problem with the modified matter potential 
V and 

/3 = exv (-4(A3i + A32)t - i 1^ dt' (V^eif) - ^F(t')) ) • (20) 
The quantity Ui is the solution of 

i^Ui = (UT?^oUo)Ui (21) 

where 

U^^oUo = ( °t ) (22) 



with 

J=^sm2ei3V^eS^ (^l^P. (23) 
A perturbative solution of Eq. ()2ip is 

Consider the integral appearing in Eq. (|24|) : 

/(T) = dtJ{t) (25) 
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= ^sin2%|^ dte-(^3i+A3.)V2exp(^-i|^*dt'(^F^.,(tO-^y(tO))^M^ 

In this integral the first exponential e~*(^3i+A32)t/2 jg rapidly oscillating function in the Sun, 
averaging to zero inside the integral. (Note that A32-R0 ~ 10^ for a 10 MeV neutrino). Indeed, 
Riemann-Lebesgue lemma ensures that 

'"''dt/(t)e^^^W ^0 



as X 00, provided that 

''dtlfit)] 



exists. The leading contribution to the integral in Eq. ()25p is 

(26) 



fit) 
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t=a 



as x — )• 00. (For a proof of this statement with integration by parts see [l3]). Using Eq. (p6|) to 
calculate the integral in Eq. (p5|) yields 



HR^) = -zsm2^,3 (A^^^.e(r = 0)) ( J ) • (27) 
For a 10 MeV neutrino with 8171^2 

~ 2 X 10^3 eV^, assuming an electron density of 100 Na / 
cm'^ in the neutrino production region, the quantity inside the parentheses in Eq. ([2] 



a = V^eir = 0) ~ 3 X 10"^ (28) 

^32 + ^31 

If one takes a = 0, one finds, e.g. on the surface of the Sun to be 

= COS^ 013^^1 + sin2 013/3, (29) 

where ipi is the electron neutrino wave function obtained in two-flavor evolution with the 
modified matter potential V. Upon averaging over the Earth-Sun distance (see e.g. [15]) one 
gets Eq. ([3]). Including the first order correction in Ui, which is proportional t a of Eq. ([28]) . 
gives 

Psxsii'e i^e) = cos^ 6li3(l-4sin^ 6*130)^2x2(1^6 i^ewithiVeCos^ 6li3) + sin^ 6li3(l + 4cos^ Oisa). 

(30) 

A different derivation of this formula was given in Ref. [16j. Even for the largest value of ^13 
quoted by T2K experiment, this correction to the survival probability is less than one percent 
for a 10 MeV neutrino. 



3. CP-violating phase 

It is rather straightforward to show that the CP-violating phase factorizes out in the neutrino 
evolution Hamiltonian of Eq. ([6]): 



Eq. (j3ip implies that the neutrino evolution operator can be similarly factorized. This 
factorization gives us interesting sum rules: Electron neutrino survival probability, -P(fe — ^ t^e) 
is independent of the value of the CP-violating phase, 6; or equivalently, in the absence of sterile 
neutrino mixing, the combination -P(f^ — ?■ Ue) + P{i't — ^ t^e) at a fixed energy is independent of 
the value of the CP-violating phase [IT]. However the probability Pii'fj, — )• I'e) depends on the 
5c p- This is why the disappearance experiments such as those that use reactor antineutrinos 
can quote a probability independent of Scp, whereas the appearance experiments, such as those 
that use accelerator neutrinos, cannot disentangle ^13 from 5cp- It is possible to derive similar 
sum rules for other amplitudes [18] . A discussion of the breakdown of this formula when sterile 
neutrinos are present is given in Ref. |19j . 

4. Concluding Remarks 

We provided a brief summary of the limits and preliminary measurements of the still poorly 
known neutrino mixing angle, ^13. Many analyses of the solar neutrino data utilize a formula 
connecting two- and three-flavor evolution. We gave a critical assessment of the validity of this 
formula. 

If the angle ^13 is rather large, there may be significant astrophysical implications. A brief 
review of the impact of the value of the mixing angle ^13 on various astrophysical phenomena 
was given in |20j . 

Acknowledgments 

I thank A. Palazzo for discussions and correcting a typo in the earlier version of the manuscript. 
This work was supported in part by the U.S. National Science Foundation Grant No. PHY- 
0855082 and in part by the University of Wisconsin Research Committee with funds granted by 
the Wisconsin Alumni Research Foundation. 

References 

[1] Nakamura K et al. [Particle Data Group] 2010 J. Phys. G 37 075021. 

[2] Balantekin A B and Yilmaz D 2008 J. Phys. G 35 075007 (Preprint arXiv7d804.3345 [hep-ph]) 

[3] Fogli G L, Lisi E, Marrone A, Palazzo A and Rotunno A M 2008 Phys. Rev. Lett. 101 141801 (Preprint 

arXiv:0806.264i [hep-ph]); Schwetz T, Tortola M A and Valle J W F 2008 2008 New J. Phys. 10 113011 

(Preprint arXiv:08O8.2O16 [hep-ph]); Gonzalez-Garcia M C, Maltoni M and Salvado J 2010 JHEP 1004 

056 (Preprint arXiv:1001.4524 [hep-ph]). 
[4] Mezzetto M and Schwetz T 2010 J. Phys. G 37 103001 (Preprint arXiv: 1003.5800 [hep-ph]). 
[5] Aharmim B et al. [SNO Collaboration] 2010 Phys. Rev. C 81 055504 (Preprint arXiv:0910.2984 [nucl-ex]). 
[6] Gando A et al. [The KamLAND Collaboration] 2011 Phys. Rev. D 83 052002 (Preprint arXiv:1009.47n] 

[hep-ex]). 

[7] Ardellier F et al. [Double Chooz Collaboration], (Preprint arXiv:hep-ex/0606025 I. 

[8] Guo Xet al. [Daya Bay Collaboration], (Preprint arXiv:hep-ex/0701029 ). 

[9] Ahn J K et al. [RENO Collaboration], (Preprint arXiv:1003.1391 [hep-ex]). 
[10] Abe K et al. [T2K Collaboration], (Preprint arXiv: 1106. 1238 [hep-ex]). 
[11] Abe K et al. [T2K Collaboration], (Preprint arXiv:1106.2822 [hep-ex]). 



H(5cp) = SU{5cP = 0)St 



(31) 



with 




(32) 



[12] Balantekin A B and Fuller G M 1999 Phys. Lett. B 471, 195 (Preprint larXiv:liep-pli/9908465l. 

[13] Balantekin A B and Yuksel H 2003 J. Phys. G 29 665 (Preprint arXiv:liep-pli/0301072 ). 

[14] Bender C M and Orszag SA 1978 Advanced Mathematical Methods for Scientists and Engineers (New York: 

McGraw-Hill) p 278 

[15] Balantekin A B 1999 Phys. Kept. 315 123 (Preprint arXiv:hep-pli/9808281 ). 

[16] Fogli G L, Lisi E and Palazzo A 2002 Phys. Rev. D 65, 073019 arXiv:hep-ph/0105080 . 

[17] Balantekin A B, Gava J and Volpe C 2008 Phys. Lett. B 662 396 (Preprint arXiv:0710.3112 [astro-ph]). 

[18] Kneller J P and McLaughlin G C 2009 Phys. Rev. D 80 053002 (Preprint .arXiv:0904.3823 [hep-ph]). 

[19] Palazzo A 2011 Phys. Rev. D 83 113013 arXiv:1105.1705 [hep-ph]]. 

[20] Balantekin A B 2010 AIP Conf. Proc. 1269, 195 (arXiv: 1006.2836 [nucl-th]). 



